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Abstract 

The thermodynamic formahsm expresses chaotic properties of dy- 
namical systems in terms of the Ruelle pressure ip{P). The inverse- 
temperature like variable /? allows one to scan the structure of the 
probability distribution in the dynamic phase space. This formalism 
is applied here to a Lorentz Lattice Gas, where a particle moving on 
a lattice of size L'^ collides with fixed scatterers placed at random 
locations. Here we give rigorous arguments that the Ruelle pressure 
in the limit of infinite systems has two branches joining with a slope 
discontinuity at /? = 1. The low and high /3-branches correspond to 
localization of trajectories on respectively the "most chaotic" (highest 
density) region, and the "most deterministic" (lowest density) region, 
i.e. ip{P) is completely controlled by rare fiuctuations in the distribu- 
tion of scatterers on the lattice, and it does not carry any information 
on the global structure of the static disorder. 

As P approaches unity from either side, a localization-delocalization 
transition leads to a state where trajectories are extended and carry 
information on transport properties. At finite L the narrow region 
around /? = 1 where the trajectories are extended scales as (lnL)~", 
where a depends on the sign of 1 — /3, if d > 1, and as (LlnL)~^ if 
d = 1. This result appears to be general for diffusive systems with 
static disorder, such as random walks in random environments or for 
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the continuous Lorentz gas. Other models of random walks on disor- 
dered lattices, showing the same phenomenon, are discussed. 

Key words : Lorentz lattice gases, chaos, thermodynamic formalism, random 
walks, localization transition 



1 Introduction 

In the past several years a large body of research has focused on the prob- 
lem of relating the macroscopic behavior of non-equilibrium systems to the 
underlying chaotic dynamics of the particles of which the system is com- 
posed. Some macroscopic transport coefficients appearing in hydrodynamic- 
like equations have been related to microscopic quantities which characterize 
the chaotic properties of the system fl], ||, ^, Ruelle, Sinai, and Bowen || 
introduced a powerful method to derive most of the interesting chaotic prop- 
erties of a given system from a free-energy type function, called the Ruelle or 
topological pressure. This thermodynamic formalism is based on a partition 
function calculated in a dynamical phase space. For systems governed by 
discrete, rather than continuous dynamics, one point in the dynamical phase 
space over t time steps consists of a trajectory, Q(t) = {xi, X2, ■ ■ ■ , Xt}, which 
is a set of t successive states of the system. The topological pressure ip{P) 
is defined as the infinite time limit of the logarithm of the partition function 
divided by the time t, in a way similar to the definition of the free-energy 
per particle in a canonical ensemble, in the thermodynamic limit. Again, in 
analogy with the methods of equilibrium statistical mechanics, there is an 
inverse temperature-like parameter j3 which allows one to scan the structure 
of the probability distribution for Q. 

This formalism has been successfully applied to Lorentz gases |]^, |^, |^, |I^] . 
These are models in which independent light particles are moving among fixed 
scatterers. They can be considered as elementary models for diffusive trans- 
port in fluids and solids. In the continuous case, the effects of disorder in the 
conflguration of scatterers can be taken into account and chaotic properties 
can be computed in the region of /5 ~ 1, as will be discussed in another paper 
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np. The model can be simplified further by constraining the light particles 
to move on a regular lattice, and placing scatterers, with some density, at 
random locations on the sites of this lattice. Such models are called Lorentz 
Lattice Gases (LLCs). For some of these we already calculated the Ru- 
elle pressure around /? = 1 in the framework of a mean-field approximation 
[§, ^, . For one-dimensional open systems on a lattice of size L we obtained 



the escape rate 7, the Lyapunov exponent A, and the Kolmogorov-Sinai (KS) 
entropy, and found good agreement with independent direct numerical esti- 
mates of the same quantities p . All these quantities depend on the average 
density of scatterers p. 

We have previously reported [|12| that, for (3 different from unity, and 
for systems with static disorder, the Ruelle pressure has unexpected proper- 
ties as L becomes large enough. In particular, in the thermodynamic limit, 
L — > 00, it becomes independent of the density of scatterers. The cause of 
this is that, for large systems, the Ruelle pressure is completely determined 
by rare localized fluctuations in the configuration of scatterers. This peculiar 
behavior is expected to be general for all diffusive systems with static dis- 
order, in any dimension. In this paper we develop the analytical arguments 
allowing for such a claim. The numerical counterpart will be presented in a 



separate paper [|I3 



In the next sections we first describe the LLG model in more detail, and 
then introduce the thermodynamic formalism. Moreover, we will extend the 
results of |]l2l to a mixed random walk model, and use it throughout the paper 
to illustrate the generality of our results. The main line of our calculation is 
to construct exact upper and lower bounds for the Ruelle pressure which, as 
we will show, coincide in the thermodynamic limit, and are determined by 
rare configurations of scatterers, except in a small region about /3 = 1. We 
refer to this situation as localization of orbits on rare fluctuations of disorder. 
This means that the dominant contributions to the Ruelle pressure in the 
limit of large systems originate from orbits (points in the dynamical phase 
space) where the particle is restricted to move on those rare fluctuations, i. e. 
for /3 < 1, on the largest compact cluster of scatterers, and for /3 > 1 on 
the largest hole. We need, as a side result, the distribution function of the 
largest cluster size over all configurations, and the crude estimate of |12] will 
be refined. The analysis of finite size effects shows that the thermodynamic 
limit is approached extremely slowly, ~ (lnL)~" where a depends on the 
model and on the sign oi 1 — (3. For finite systems, we have estimated the 
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/9-rangc around /9 = 1 in which the Ruelle pressure is still determined by 
trajectories extending over the whole system. As (3 is deviating more and 
more from unity, the orbits become more and more localized on the largest 
cluster or in the largest hole of the entire configuration. In one dimension, 
there is an intermediate state with "weak" localization ( see section 7). 
The extension to continuous Lorentz gases is presented in section 8. 



2 Lorentz Lattice Gases 

A 'light' particle moves ballistically in a finite simple cubic domain V hav- 
ing periodic or absorbing boundaries and containing V — L'^ sites of a d- 
dimensional cubic lattice. The allowed states of the system x = {r, Cj} at 
time t (t = 0, 1, 2, • ■ ■) are specified by the position r G P and the velocity 
Cj of the moving particle. The set of possible velocities Cj {i = 1, 2, ■ • ■ , 6) 
connects each site to its b nearest neighbors, where the coordination num- 
ber b equals 2d for a simple hypercubic lattice. A fraction p of the sites 
- chosen at random - is occupied by a scattcrcr or 'heavy' particle. The 
quenched configuration of scatterers is specified by the set of Boolean vari- 
ables {p(r), r e V}, where p(r) = 1 if site r is occupied by a scatterer, and 
p(r) = if site r is empty. When the hght particle hits a scatterer, it is 
scattered to one of the lattice directions with a probability that depends on 
its incident velocity. The scattering laws arc further specified by introduc- 
ing a transmission coefficient p, a reficction coefficient q, and, for hypercubic 
lattices, a defiection coefficient s, normalized as 

p + q + 2(d-l)s^l. (1) 

More formally, Wij with i, j = {1, 2, • • • , 6} is the probability that the moving 
particle with incident velocity Cj is scattered to the outgoing velocity Cj with 
normalization J2i = 1. For instance, on a square lattice, the transition 
matrix has the form 
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The scattering at site r is described by the random transition matrix 
Wij{r) which depends on the configuration of scatterers {p(r) ; r e V}, and 
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is given by 

m,{v)=p{v)W,, + {l-p{v))5,,. (3) 

At full coverage (p = 1) the moving particle performs a random walk with 
correlated jumps, referred to as the Persistent Random Walk (PRW) |T^. 

The time evolution of this system, in a fixed configuration of scatterers, is 
described by the Chapman-Kolmogorov equation for the probability 7r(x, t), 
with X = {r, Cj}, to find the moving particle at time t on site r with incident 
velocity Cj, i.e. 

7i{x,t + l) = J2w{x\y)n{y,t). (4) 

y 

In the case of absorbing boundary conditions, boundary states y = {r',Cj} 
referring to a particle entering the domain V are excluded from the y- 
summation. The transition matrix w{x\y) represents the probability to go 
from state y = {r', Cj} to state x = {r, Cj}, and is given by 

w{x\y) = 6{r-c,,r')W,,{r'). (5) 

The basic ideas of this paper are applicable to the much wider class of 
diffusive models with static disorder, such as hopping models with bond or 



site disorder [0, as well as to continuous Lorentz gases (see section 8). As 
the most immediate generalization of a LLG, we consider another model of a 
random walk, called a mixed random walk (MRW), in which a particle moves 
on a lattice filled by a random mixture of X types of scatterers. This model 
may be described by X scattering matrices of the form of Eq.(^, i.e. W^^ 
with parameters Pk,(lk,Sk ^ k < X). The model contains the 'ballistic' 
LLG, described above, as the special case with X = 2, pi = 1 — qi = p, and 
P2 = 1 — '?2 = 1 or wj;p = 6ij. The scattering at site r in the MRW- model is 
then described by the random transition matrix 

W,A^) = J: Pk{r)wl'\ (6) 

k=l 

where /)fc(r) = 1 if site r is occupied by a scatterer of type k, and zero 
otherwise. 

Boundary conditions may be either periodic (closed system) or absorbing 
(open system) on the boundaries of domain V, and the transition matrix 
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satisfies the normalization relations 

The inequality sign in Eq.(|^) for open systems refers to the case where 
y = {r, Cj} denotes a state at a boundary site r with non-entering veloc- 
ity (boundary states with entering velocity do not occur). Indeed, the sum 
over X excludes states where the particle has escaped from the domain V. 
Hence the probability for remaining inside the domain decreases when the 
particle finds itself on a boundary site. 

3 Thermodynamic formalism 

As stated in the introduction, the starting point for this paper is a partition 
function defined in the dynamic phase space whose points Q{t) represent 
trajectories of t time steps : 

^L(/3,t|xo)=E[^(^'^l^o)]^ (8) 
n 

where P{Q,t\xo) is the probability that the system follows a trajectory 
Q{t) = {xi, X2, ■ ■ ■ ,Xt}, starting from xq at t = in a given system of linear 
dimension L = V^l'^. The temperature-like parameter /3 allows us to scan the 
structure of the probability distribution P, where large positive and negative 
/5-values select respectively the most probable and most improbable trajecto- 
ries [H. The concepts used in this section have been discussed in great detail 
in Ref.ia. 

In each specific system the probability P(f2, t|a;o) for a given trajectory 
can be expressed in terms of the transition probabilities w[x\y): 

P(n,t|a;o) = n^=iti;(x„|x„_i). (9) 

The partition function is determined by the properties of the matrix wp[x\y) = 
[w{x\y)]^, which is defined by raising each matrix element w{x\y) to the power 
f3. For large times the partition function for almost all systems becomes in- 
dependent of the initial point xq (ergodicity, see JlOl), and is determined by 
the largest positive eigenvalue Al{i3) of the matrix Wi3{x\y), which for ergodic 
systems can be shown to be non-degenerate. 
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There is a slight comphcation because hypercubic lattices are bipartite, 
i.e. the moving particle is always on even sites at even times and on odd 
sites at odd times, or vice versa. The system therefore consists of two inde- 
pendent ergodic components, which should be considered separately, and the 
matrix Wf3{x\y) is called a periodic matrix of period two To avoid these 
complications one may consider the time t to be an even integer multiple of 
the time step and then consider the matrix w^{x\y) = J2z'^i3{^\^)'^/3{^\y)j 
defined between even or between odd sites only, to be the fundamental ma- 
trix. Again, the largest eigenvalue [Ai(/3)]^ of the matrix Wp{x\y), restricted 
to one sublattice, is non-degenerate. 

In addition, the topological or Ruelle pressure is defined as 

(/?, p) = lim I (In Zl {[3, t\xo))p , (10) 

where (. . .)p denotes an average over all configurations generated by the pre- 
scription that for each lattice site, independently, p is the probability that 
it will be occupied by a scatterer. The topological pressure is independent 
of xo, if the system is ergodic, and can be expressed in terms of the largest 
eigenvalue, Al(/3), of the matrix wp as 

^i(/5,p) = (ln(Ai(/3)))p, (11) 

where we have taken the infinite time limit inside the configurational average. 

Several chaotic quantities can be derived from this function [||, 0. For 
example: the sum of all positive Lyapunov exponents is A = Ya^'^ \ = 
—ip'i^^l); the escape rate for open systems is 7 = —ipL^l); the Kolmogorov- 
Sinai entropy follows from the generalization of Pesin's theorem to h^s = 
— V'l(I); the topological entropy hx satisfies hx = ^l(O); the Hausdorff 
dimension dn of the repeller (the set of trajectories that never escape) for 
an open system is the zero-point of the Ruelle pressure, i.e. ipLidn) = 0. A 
prime in the above formulas denotes a /3-derivative. 

4 Upper and Lower Bounds 

In this section the Ruelle pressure will be calculated in the limit of infinite 
system size, by constructing upper and lower bounds at finite L and analyzing 
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their limiting behavior for large L. Consider first the Lorentz lattice gas. For 
a closed system it follows from the definition of Wf3{x\y) and Eq. (H) that 

Y: wpix\y) = pir)Wi^) + (1 - p(r)) (12) 

X 

with y = {r, Cj} and 

W{(3) =a + b + 2{d-l)c 

a=p^, b = q^- c=s^ (13) 

For open systems the equality sign in Eq.(|l^) is replaced by a 'less than' sign 
in case y is a. boundary state. As a general upper bound, valid for both open 
and closed systems, we obtain 

'^Wfs^xly) <W{P) (r = scattering site) 

X 

'^wp^xly) <1 (r = empty site). (14) 

X 

If /3 < 1, this implies that Y^x'^pi^lu) < everywhere, as W{P) > 1. 

This inequality combined with Eqs.(||) and (^ yields: 

ZL{(3,t\xo) < W{[3)ZL{f3,t-l\xo) < {W{f3)y. (15) 

Then the pressure, defined by Eq.([To|), satisfies the inequality: 

MP,P)<^^W{P). (16) 

If P > 1, and consequently W{P) < 1, the analog of Eq.(|^) becomes : 

ZLi/3,t\xo)<l (17) 

and 

The above equations, (16) and (18), provide upper bounds for the Ruelle 
pressure for all /3-values. 

To construct a lower bound to Zl we consider clusters of scatterers, i.e. 
regions where every site is occupied by a scatterer, and select the cluster of 
largest size M. A cluster is said to have size M, if the largest inscribed cube. 



8 



oriented along the lattice directions, has a linear dimension M. If there are 
several largest clusters of the same size, just choose one arbitrarily. The value 
of M is well defined for any given configuration of scatterers. 

As all terms in the sum (H) are non-negative, any sum over a subset of 
trajectories will give an exact lower bound for Z^. Let Z^(/3,t|xo) denote 
the sum over all trajectories which remain confined for t time steps to the 
largest inscribed cube of size M (again, if for a given cluster there is more 
than one inscribed cube of linear size M, choose one arbitrarily), then we 
have a lower bound: 

zl > ^r- (19) 

In fact, ZfJ^ is equal to the partition function of a PRW in an open hyper- 
cubic domain with M'^ sites. According to Eq.([TT|) this requires the largest 
eigenvalue Af)^(/?) of the matrix Wi3{x\y) for the PRW, which can be found 
in Ref . |]1 7|| , and reads for sufficiently large M- values : 

Ar(/5) = - A{f3)e + Oik')} (20) 

where fc^ = X]a=i 

Here k is the smallest wave number accessible to the system, i.e. A; = for a 
closed system (with periodic boundaries) and k^ — n/M {a = x,y, ■ ■ ■ ,d) for 
an open hypercube (with absorbing boundaries). On the basis of Eqs. (p!l|) 
and (|T^) we find the first lower bound on the Ruelle pressure, valid for all 
/?- values: 



M(3,P) > (lnAr(/5)) 



p 



= \nW{/3)- A{/3){dTT'^/M^)p. (22) 

If we can show that the moment < > tends to zero when the system 

size increases, then this lower bound will tend to the upper bound (|TB|) in 
the range (3 < 1. In the range (3 > 1, this lower bound will approach the 
finite negative value hiW{(3). 

We need another lower bound which will tend to the upper bound (0) for 
/3 > 1. In order to find it, we consider for any fixed configuration the longest 
line segment free of scatterers. Contrary to the largest cluster defined above. 
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the largest empty line segment is always a one-dimensional domain, whatever 
the dimensionality of the system is. Let M be the number of empty sites 
on this line segment, and Z the partition sum (H) restricted to trajectories 
confined to this line segment. In fact, we keep only a single trajectory, 
which travels continually through the empty region and is reflected by the 
two scatterers at the end sites. For sufficiently large times the number of 
reflections is approximately t/M. 

The resulting sum is Z ^ g^*/*^ and we have thus a second lower bound 
for the Ruelle pressure, valid for all /3-values, 

V'L(/3,p)>/?(lng)(l/M),. (23) 

In summary, the following upper and lower bounds apply to all LLCs : 

\TiW{(3)-A{(3){dTx^/M'')p < M(^.P) <\nW{(3) (/?<!) 

(5{hiq) (1/M), < MP.P) < (/?>!) (24) 

The above bounds can be extended straightforwardly to the mixed random 
walk (MRW) models, where Eq. (|12D becomes 

Y.wMy) = j:h{^)w^'\(i)- (25) 
X k=i 

with W^''\P), ak,bk, and Cfc defined in a similar way as in Eq.([TB|). 

The lower bounds and Z^ are respectively determined by the 

largest cube of PF"'"-scatterers containing M'^ scatterers, where is the 
type of scatterers for which InW^'^^P) is largest for a given /3. The resulting 
upper and lower bounds in MRW-models can then be summarized as: 

lnW+{(3) - A+{f3){d7TyM^)p < ^l{P,p) < \nW+{(3). (26) 

The bounds for /3 < 1 contain the LLG as a special case; the bounds for 
/3 > 1 are different. 

5 Thermodynamic Limit 

The goal of this section is to show that the upper bounds of section § are 
indeed the asymptotic values for the Ruelle pressure in the limit of infinite 
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systems {L oo). To do so, we need to evaluate the inverse moments, 
< M~*^ > (A; = 1,2), entering in Eq.(p^, in the hmit as L — cx). This 
requires the asymptotic behavior of the probabihty that the largest cluster 
is of size M. 

We first consider the one- dimensional case where configurations are gen- 
erated by distributing scatterers on the lattice sites according to the prescrip- 
tion that the occupation probability for each lattice site is p, independently of 
the other sites. Then the total number of scatterers N may fluctuate around 
its average value pL. A crude estimate can be obtained by noticing that the 
average number of clusters of size m is approximately Lp™(l — pY ■ Indeed 
the cluster can be centered on L different positions on the lattice (or L — m 
positions for an open system), it contains m scatterers, and is bordered by 
two empty sites. For m to be a typical value for the size of the largest cluster, 
the above expression must be of order unity, which implies that M scales as 
InL. For the inverse moments of M this implies < M"'^ >~ (InL)"'^ for 
large L. Hence upper and lower bounds in Eqs.(^) approach the same limit. 

This argument can be extended directly to higher dimensions. A cluster of 
size m (this means that the largest inscribed cube has side m) occurs roughly 
l^dpTn i^jj-Qgg^ where we used that for large m the probability of finding at least 
one empty site in each of the boundary hyperplanes is very close to unity. 
For L sufficiently large LV'" is of order unity if ~ InL. Consequently 
the inverse moments < M"*^ > ~ (InL)"'^/'^ for L oo. 

The L-dependence of the inverse moments can be obtained more rigor- 
ously by the following observation. We identify the clusters of size m, with 
m ^ 1, as 'noninteracting molecules' of species m with partial densities 
{n{m) ~ ; m = 1, 2, ...}. The probability to find a volume of size V = L'^ 
unoccupied by clusters of size > M is then 



P(M) = exp 



exp 



~V ^ n(m) 

m>M 

-L'p'n . (27) 



Here we replaced the sum in the exponential by the first term, since the 
size of subsequent terms in the series decreases extremely rapidly. We also 
replaced M + 1 by M, which will induce some correction terms of relative 
order (InL)"^/'' in the final expression (|30|) . The probability that the largest 
cluster is exactly of size M is then A{M) = P{M) - P{M - 1), or in the 
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continuum limit, 



A{M) = P'{M) ~ -dL'^ ln(p) M'^-^p^^ P{M). (28) 

For large L the inverse moments < >= dMM~^A{M) can be eval- 
uated asymptotically by a saddle point method, as A{M) is sharply peaked 
around its maximum. The maximum is located at Mq, which is the root of: 

[In A{M)]' = -rfL<^(lnp)M"'-V*^' + c/ln(p)M'^"i + {d-l)M-\ (29) 

For large M the solution of this equation is determined by the first two terms 
on the right hand side (dominant balance argument), yielding, 

with correction terms of relative order (lnL)~^^'^. For large L the inverse 
moments behave asymptotically as, 

<M-'>~M-=(^)""' (31) 

in agreement with the crude estimate above. 

In the preceding paragraphs the probability p of occupation of sites by 
scatterers has been kept fixed. It can be shown that the results (^0]) and 
( PID for the L dependence are still correct if one fixes the total number of 
scatterers, N, in all configurations 



The above results can also be used to estimate the typical size M of 
the largest empty line segment. The distributions for scatterers and holes 
are symmetric by exchange of p and 1 — p. As we are interested in a one- 
dimensional domain, whatever the dimension d of the lattice, we have to 
replace M'^ by M. We straightforwardly obtain 

<l/M>^l/M„.-!^. (32) 

We conclude therefore that in the thermodynamic limit L ^ oo, the moments 
< M"*^ > vanish and thus the lower bounds of section ^ converge towards the 
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-0.2 I ■ ' ■ ' ■ ' ^= 

0.0 0.5 1.0 1.5 

parameter p 

Figure 1: Ruelle pressure for a MRW in a one- dimensional system of infinite 
size for pi = 0.9 and p2 = 0.7. It is independent of the density of scatterers 
as long as pi or p2 is not or 1. Note the difference with the LLG, where 
the branch for (3 > 1 is = 0. 



L-independent upper bounds. This yields for the Ruelle pressure in LLCs 
in the thermodynamic limit 

2im^,(/3,p) = | Q ^^^^^^ (33) 

In the MRW-models, the Ruelle pressure is 

\im iIjl{(3,p) = In W+ {(3) (34) 

L— >oo 

as is illustrated in figure ^ for MRW models. 



6 Localization, extension to MRW's 

In the previous section we have shown that the dynamic partition function 
and the Ruelle pressure of LLCs in thermodynamically large systems are 
completely determined by the rare fiuctuations in the spatial distribution of 
scatterers. It is worth stressing again here that for the regions (3 < 1, (3 = 1, 
and (3 > 1, and L oo, different sets of trajectories make the dominant 
contributions to the Ruelle pressure. For (3 < 1 the trajectories contributing 
to the partition function in a given configuration of scatterers are localized on 
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the largest compact cluster of scatterers, i.e. localized in the 'most chaotic' or 
'least deterministic' region. For f3 > 1 only a single trajectory contributes, 
which — on any (i-dimensional lattice — is localized on the longest line 
segment that is free of scatterers, i.e. the relevant trajectory is localized in 
the 'least chaotic' or 'most deterministic' region of configuration space 

Therefore, as L — ^ oo the Ruelle pressure becomes independent of the con- 
figuration (except for atypical configurations such as strictly periodic ones); 
it is even independent of the density of scatterers (except at p = 0, where 
no fluctuations exist any more). It does not carry any information on the 
structure of the random medium. 

Since for (3 different from unity, relevant trajectories do not explore the 
whole system but only a small part of it, the 'mean field configuration' with 
all scatterers more or less equidistant from one another is not at all a typical 
configuration. On the contrary, it is the one that gives, among all configu- 
rations, the minimal value for the Ruelle pressure. On the other hand, the 
maximal Ruelle pressure is obtained for the configuration where all of the 
scatterers form a single compact cluster. The average will be somewhere in 
between. This means that any calculation starting from a mean field approx- 



imation will give very poor results for /5-values different from unity |23 



However at /? = 1, the Ruelle pressure for the LLG and its derivatives 
with respect to (3 do depend on the overall density M and on more details 



of the total configuration of scatterers [ 15 1. For finite L, this is also true in 



a small region around P = 1. There, relevant trajectories are extended or 
delocalized, and explore large regions of the lattice. This conclusion is based 
on the reasonably good agreement for escape rates and Lyapunov exponents 
between the results from computer simulations and mean field calculations 
for the LLG §. 

The same conclusions carry over to the MRW-models, where in the ther- 
modynamic limit the trajectories are localized on the largest compact cluster 
with iy"'"(/5)-scatterers. At (3 = 1, all scatterers have the same InW = 
and again trajectories are delocalized on the whole lattice. The structure of 
the typical mean field configurations, contributing around (3 = 1, has not yet 
been explored, and mean field estimates for the Lyapunov exponents in open 
systems have not yet been derived for MRW's. 

Interesting new phenomena can occur near those values of /3 7^ 1 where 
different types of scatterers may have the same value of In W. For example, 
this occurs at = if all scatterers have the same number of non-zero- 
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scattering directions. Then the moving particle cannot distinguish between 
the different types of scatterers and the relevant trajectories become again 
"delocalized" on a large cluster with a random mixture of the different scat- 
terers with the same value of lnW{(3). More explicitly, if there are X types 
of scatterers, it may occur that K {2 < K <X) of them have the same In 
strictly greater than the In W for all other types of scatterers. Then the clus- 
ter of ly^-scatterers on which the relevant trajectories are localized contains 



a random mixture of these K types of scatterers |15 



Suppose now that we consider, for (i > 1, a LLG or a MRW model that 
shows a percolation transition. In such a case is important to note that 
the definition for the cluster size used here is not the number of connected 
sites, but the size of the largest inscribed hypercube, which typically is much 
smaller than the system size even for a percolating cluster. Thus the per- 
colation transitions in such models have no effect on our considerations and 
the results of this paper remain valid. 



7 The delocalization region 

In this section we estimate the size of the delocalization region around (3 = 1 
in LLCs for finite systems. The Hausdorff dimension of the repeller (i.e. the 
set of trajectories which do not escape from the system after an infinite time) 
is a root of the Ruelle pressure, 

^ddn) = 0. (35) 

For a large system, dn is close to unity 0. Using the facts that il^hiX) = 
—7 ~ Ddir'^/LP' for a hypercubic domain in (i-dimensions, where D is the 
diffusion coefficient, and V'l(I) = ~H\^>q\ — ~^oq where Aoo represents 
the sum of positive Lyapunov exponents in the infinite L-limit, we find that 
the Hausdorff dimension dn for a large hypercubic domain is in first approx- 
imation, 

dHc^l- I — ]—, (36) 



where D and Aqo depend on the density of scatterers. Therefore, as the 
structure of the repeller is a fundamental feature of the system, the crossover 
region should extend at least over a /3-range of order On the other hand 
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we have concluded that for L ^ oo and for (3 different from unity, the Ruelle 
pressure becomes independent of the global structure of the disorder, as 
the relevant trajectories become localized in regions of the lattice where rare 
fluctuations of high or low density of scatterers occur. This was demonstrated 
in previous sections in the limit of infinite systems. Therefore as long as 
the mean field value of the Ruelle pressure or of the largest eigenvalue is 
smaller than the lower bound, the states of the system are localized. In fact, 
numerical results [1^ support our intuition that the effect of localization can 



be estimated fairly well for any large but finite system by taking the lower 
bounds in Eqs. (p^ as estimates of the Ruelle pressure for /? < 1 or > 1. 
Hence, for a given L crossover from localized to extended states occurs as we 
approach (3 = 1 from either side, when the mean field value equals the lower 
bound. 

To obtain an estimate of these crossover values, we compare the Ruelle 
pressure of the mean field configuration with the estimate for the Ruelle 
pressure based on the lower bounds (see Eq. (p2D ). It is equivalent to compare 
the eigenvalues of the matrix wp associated with a localized and with a 
delocalized eigenstate. 

The second one is obtained for the 'mean-field' configuration in LLCs 
from the PRW expression (|20|) by a rescaling argument H, |13|. It reads 

Ar(/?) = {Wm' (l - pA(/3)d(7r/pL)2) + O^l/L') (37) 



with W{I3) and A defined in Eqs.(|T3D and (|2T|). 

For (3 < 1 localization takes certainly place if the lower bound on the 
Ruelle pressure is larger than the mean field value, i.e., 

hiW{f3)- /\{l3)dTi^/Ml > ln{W{p)Y + 0{1/L^), (38) 

where the left hand side is the lower bound given in Eq. ( P^ with < >~ 
Mg"^ on account of (0). The right hand side is the mean field value given 
by Eq. (PTj). By expanding both sides in powers of e = 1 — /3, we find that 
the Ruelle pressure is determined by localized trajectories only if 



e > e_ 



6{l-p)Mi S{l-p) 



d\nL 



2/d 

(39) 



where 6 = \p\np + glng + 2{d — l)slns| > 0. We note that c?// = 1 — e in 
Eq.(|36|) is indeed within the delocalized region, as was to be expected. 
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However, the crossover between localized and delocalized states may in- 
volve some intermediate states. In principle one might have a "weak localiza- 
tion" in a region of size Mi (with In L <^ Mi <^ L) , where the local density, 
p + Ap, is slightly larger than p, but where Ap is large enough so that tra- 
jectories, remaining confined to this region, dominate the Ruelle pressure. 
To estimate the largest density fiuctuation to be found in a region of size 
Ml we first note that the probability for a density fiuctuation Ap in such 
a region can be estimated as exp(//ApMf), with p the chemical potential 
of the scatterers (considered as lattice gas particles). Since the number of 
different regions of this size is on the order of L'^ for Mi in the above range, 
the largest density fiuctuation occurring in one of these regions follows from 
the requirement L exp or 

-pApMf~lnL^. (40) 

In a manner similar to ( ^Sf ) we compare the mean field value ~ In {W{l3)Y 
with the mean field estimate of the Ruelle pressure corresponding to trajec- 
tories confined in a region of average density p + Ap, 



^ff^ ^ In [W{(3)r^'' 



A(/3) 



(p + Ap)M2 



(41) 



Expressing Ap in terms of Mi according to (|40D and expanding the difference 
'^/'Aff ~ ^l'^ lowest order in e yields a condition for "weak localization" in 
the most dense region of size Mi, namely. 



^l^f - " - - > 0. (42) 



■•MF „i.MF ^ 5d\nL _ (iA(l)7 r^ 
~^ Mf pMl 

By taking the estimate of the delocalization region in ( pOf ) one immediately 
sees that for e ^ (InL)^^/'^, this inequality can only be satisfied in c? = 1. 
Localization will occur on a region of size Mi maximizing the difference (H3) 



in Ruelle pressures. Differentiating (|42|) with respect to Mi gives an Mi 



that is proportional to l/(elnL). As long as this is -C weak localization 
will occur. As soon as e ~ = 2A7r^/(p(5 L In L) the confinement region of 
the dominant trajectories becomes comparable to the full system and weak 
localization is no longer a meaningful concept. Hence, in one dimension one 
can distinguish in addition the weak localization regime e^^, < e < e_. The 
region of full delocalization is narrowed down to 1 — (3 < eyj. 



17 



For /3 > 1 and d > 1, the crossover value can be determined by comparing 
the mean field estimate with the lower bound (|23|) for the Ruelle pressure in 
the LLG, combined with Eq. (PT|), i.e. 



For d = 1 one can show again, by using arguments similar to those above, that 
there is a region of weak localization where the density is slightly lower than 
average. This region occurs for /3 values given by C"/(L In L) < P — 1 < e'_^_, 
where C is a positive constant. 

In summary, the Ruelle pressure in LLG models is determined by ex- 
tended or delocalized states if P is in the interval {1 — a_/(ln L)^/*^, 1 + 



a+Zilniy/'^} , for d > 1, and in the interval {l-C/{L\nL),l + C' /(LlnL)} 



for d = 1, where a±, C, and C are some positive constants. 

8 Extension to continuous systems 

Our considerations can be extended to the case of continuous Lorentz gases 
with static disorder. For jS < 1, we expect the moving particle to be localized 
in a region of space with a high density of scatterers, while for /5 > 1 it 
should be localized in a large region where the density of scatterers is zero. 
To provide some qualitative explanations of this observation, we consider 
a Lorentz gas with hard spherical scatterers, the so-called non overlapping 
Lorentz gas. Extension to overlapping scatterers or soft scatterers is possible 
but will not be considered here. In the non-overlapping Lorentz gas and 
(3 < 1 the Ruelle pressure approaches in the thermodynamic limit that of 
a closely packed system of hard spheres of diameter a. To understand this 
it suffices to bound the dynamical partition function by the contribution of 
all trajectories confined to a hypercubic volume of size M"^, containing the 
centers of Af = pM'^ scatterers. The probability of finding such a volume in 
the system can be estimated conservatively to be at least proportional to 




(43) 



By expanding in powers of /5 — 1 = e' we find localization for 




(44) 




(45) 
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By expanding the logarithm of the ratio of configurational partition functions 
in powers of the volume M'^ of the hypercube (using M = pM'^) — as the 
hypercube is a small subsystem of the total system with (A^, V) — and noting 
on the other hand that the partition function Q{^/, M'^) of the hypercube 
increases exponentially in M'^, we conclude that the probability in ( ^5|) is 
proportional to {V/a'^) exp(— aikf^) with a some constant. 

For any average density below the close-packing density this can be made 
of order unity by choosing M'^ proportional to ln(V^/a'^), which implies that 
for increasing V arbitrarily large volumes with a density arbitrarily close to 
the close-packing density can be found. 

It is not clear what happens to the Ruelle pressure when all of the particles 
can move, as in fluids, for example, although it seems obvious that for (3 > 
1, trajectories where the particles rarely collide will dominate the Ruelle 
pressure. 

9 Conclusion 

We conclude this paper with a number of remarks: 

1) In this paper we have discussed the Ruelle pressure for diffusive models 
with static disorder. Our results indicate that for large systems and for all 
but a small range of values of the inverse temperature-like parameter j3, the 
Ruelle pressure is determined by rare fluctuations in the configuration of 
scatterers, and consequently carries no physical information on the chaotic 
scattering of the moving particle during its motion through the frozen-in 
disorder. Only in a narrow region around /? = 1 does the thermodynamic 
formalism yield physically relevant information on the chaotic scattering in 
diffusive systems with static disorder. The extension to continuous systems 
outlined in has been made explicit here without having to use the usual 
tools of kinetic theory. The localization phenomena in the Ruelle pressure are 
in fact most similar to the asymptotic behavior (stretched exponential decay) 
of the survival probability of a random walk in a random array of absorbing 
traps. The survival probability is solely determined by the extremely rare 
fluctuation that the random walk flnds itself in the largest region free of traps 

In different areas of statistical physics analogous phenomena occur, where 
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the large time or the small frequency /energy asymptotics are controlled by 
extremely rare spatial fluctuations, such as in Lifshitz tails [20], Griffith's 
singularities and directed polymers . 



2) It has been shown numerically in one dimension |T3[ that for finite sys- 
tems and outside the crossover region the lower bound found for the Ruelle 
pressure is also a good estimate for the pressure itself, indicating that local- 
ization on the largest cluster indeed occurs. We conjecture that in higher 
dimensions it is the largest convex cluster inscribed in a set of connected 
scattering sites that will determine the Ruelle pressure. 

3) To allow the dynamical partition function to scan the full structure 
of our diffusive models with static disorder, time should be sufficiently large 
that the moving particle can explore the entire volume of the system. Conse- 
quently t » L^, which determines the physically relevant order of limits. In 
determining the Ruelle pressure in Eq. one takes first the limit t ^ oo for 
fixed L, and next allows L to tend to oo. Therefore, for a fixed system size 
L, the trajectory has an infinite time to explore the system and to find the 
largest cluster, where it will then stay localized with a high probability. An 
interesting open problem remains to explore both the time and size depen- 
dence of the dynamical partition function, Eq. (8), to see how the various 
features discussed here are approached in the limit of infinite time, but finite 
size; to study diffusive behavior when the time is kept finite and the size of 
the lattice is allowed to become infinite; and to study the behavior of the 
dynamical partition function when both L and t approach infinity in some 
coupled manner. 
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